L-bl  au 


TECHNICAL 
MEMORANDUM 


No.  K-88/63 

VARIATION  OF  PARAMETERS 

BY  POISSON’S  METHOD 
G.  E.  0.  Giacaglia 
Computation  and  Analysis  Laboratory 


D  D  C 


!7T'r’ 


FEB  88 TO 

U  L 


D 


U.  S.  NAVAL  WEAPONS  LABORATORY 
DAHLGREN,  VIRGINIA 


U.  S.  NAVAL  WEAPONS  LABORATORY 


TECHNICAL  MEMORANDUM 
December  1963 


No.  K-88/63 

VARIATION  OF  PARAMETERS 
BY  POISSON* S  METHOD 
G.  E.  0.  Giacaglia 
Computation  and  Analysis  Laboratory 


Approved  by: 

RALPH  A.  NIEMANN 
Director,  Computation  and 
Analysis  Laboratory 

This  memorandum  is  not  to  be  construed  as  expressing  the  opinion 
of  the  Naval  Weapons  Laboratory,  and  while  its  contents  are  considered 
correct,  they  are  subject  to  modification  upon  further  study. 


TABLE  OF  CONTENTS 


NWL  Technical  Memorandum  No*  K-88/63 
Table  of  Contents 

Page 

ABSTRACT  1 

FOREWORD  1 

INTRODUCTION  1 

Integrals  of  Motion  1 

POISSON1 S  VARIATIONAL  FORMULA  2 

Application  to  Celestial  Mechanics  3 

APPLICATIONS  3 

Integral  of  Energy  3 

Angular  Momentum  4 

Integral  of  Perlcenter  5 

Time  of  Perlcenter  8 

SUMMARY  8 

REFERENCES  9 

****** 

Distribution 
Mj  (Code  5130)  (2) 

BUWEPS  (HT  Division)  (2) 

V 

K-l 

KA  (20) 

KP  (2) 

ACL  (5) 


ABSTRACT 


A  direct  derivation  for  the  variation  of  the  elliptic  elements  is 
given,  starting  from  the  variational  equation  due  to  Poisson.  The 
method  is  suitable  for  any  integral  from  elliptic  motion  and  gives 
the  variation  of  an  element  by  the  use  of  a  single  equation. 

FOREWORD 


This  work  was  undertaken  in  support  of  NRL  Work  Request  No.  SPASUR 
WR- 3-005.  The  purpose  of  this  report  is  to  give  a  simple  method 
for  deriving  the  equations  of  variation  of  constants  associated 
with  Keplerian  motion. 

INTRODUCTION 

The  well  known  method  of  variation  of  arbitrary  constants  for  the 
solution  of  a  linear  system  of  differential  equations  is  usually 
applied  in  celestial  mechanics  by  choosing  as  the  arbitrary  constants 
a  set  of  elements  of  the  Keplerian  unperturbed  motion.  The  present 
memorandum  shows  how  the  same  method  can  be  used  by  taking  the 
arbitrary  constants  equal  to  the  first  integrals  of  the  unperturbed 
Keplerian  motion.  In  some  cases  this  leads  to  a  simplification  of 
the  formulae  needed. 

Integrals  of  Motion 

We  will  consider  first  elliptic  unperturbed  motion,  though  the 
method  is  suitable  for  any  kind  of  motion  under  a  central  Newtonian 
attraction. 

1*  Energy 


2.  Angular  momentum 

-*  i  -»  -• 

C  =|r  x  v 

3.  Integral  of  perlcenter  ( Laplace) 

-4  -4 

p  =  ue  i 

where  i  is  a  unit  vector  pointing  toward  pericenter,  from  the  center 
of  attraction. 
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4.  Time  of  pericenter 


where  /is  the  mean  anomaly. 
POISSON’S  VARIATIONAL  FORMULA 


Consider  a  system  of  differential  equations  reduced  to  its  canonical 
form: 


x^  =  fj^Xjt)  (i  —  1,  2,  ««.,  n)  (l) 

and  let  X(x,t)  =  C0  be  a  first  integral  of  this  system,  called 
"unperturbed."  We  are  using  x  for  the  set  (x1#  x2,  . ..,  xn). 

Consider  now  the  "perturbed"  system: 

xi  “  fi (x,t)  +  8i(x,t)  (i  =  1,  2,  ...,  n)  (2) 

The  problem  is  to  find  the  variation  of  C0  due  to  the  addition 
of  gj(x,t).  Suppose,  therefore 

X(x,t)  =  C(t) 


From  this 


c  =  v  x, 

at  j  aXj  J 

=  —  +  £  [f  j(x,t)  +  gj(x,t) 
St  JSxjL 


£  —  g,(x,t) 

j  aXj  J 


Therefore 


C  =  E  g,(x,t) 
j  dXj  J 


(3) 
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Application  to  Celestial  Mechanics 

The  formula  (3)  is  now  applied  to  the  system 


“4  —4  1 1  /  \ 

r  =  v,  v  =  -  iij.  r  +  R  (4) 

where  ^  is  the  disturbing  force  per  unit  mass.  This  system 
corresponds  to  six  scalar  differential  equations  in  the  space 
(x,  y,  z,  vx,  vy,  vz).  We  know  first  integrals  when  ^  =  0. 

In  the  case  where  equation  (4)  applies. 


C  =  <p(r, 


and 

c  =  ^E-Rx+^_Rv+^2-R,  =  grad-  <p  .  R, 

3vx  Svy 

so  that 

fi-R  grad^  <p  (Poisson).  (5) 

Relation  (5)  is  quite  general,  C  being  any  integral  of  the 
unperturbed  motion. 

APPLICATIONS 

1.  Integral  of  energy 

h  =  I  v!  -  t 

2  r 


h  =  gracL. 

v 


R  =  v 


Now,  in  a  Keplerian  motion 


h  = 


-  li- 
Sa 


h  -•  JL.  a 
2a2 
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and  therefore: 


*  Oq2  -*  -♦ 

a  =  —  i  v  •  R 
4 

2.  Angular  momentum 

C  =  r  x  v  =  Ck 

with  respect  to  an  equatorial  system: 


(6) 


sin  ft  sin  I 
k  =  [  -cos  fi  sin  I 
cos  I 


Therefore 


Note  that 


5  =  dT  +  c(|n  +  |i) 


k  •  k  =  1 


m- 


sin2  I 


Ic  •  =  0 


£.S.o 

dn  di 


=  o 

bl 


Thus  we  obtain 


«4  • 

C  •  k  a  c 


t*  V  =  X 
bl. 


m  «•* 

C'^sC  sin2  I  h 

b(l 


c.|ksci 

ol 
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Applying  (5) 


and  then 


C  =  r  x  R 

•  -♦  -4  — ♦ 

C  =  (r  x  R)  •  k  =  tr,  R,  k] 


If  we  consider 


then 


2ae  jiae 


Using  (6)  and  (7),  there  results 


Let  us  deduce  p  as  function  of  r,  v.  Consider  the  quantity 


Prom 


v  =  -  14—  r .  r  x  v  =  C 
r3 


there  results 


c  X  7  =  Ugfx  (?  x7)  =  [(?  .  7)  ?  -  ?j 

=  „  [r  .1- r  v  1  =  -A  (a  l  V  _i  (C  X  v) 

L  r2  J  dt  \  r  y  dt 

Therefore 

?  x  v  =  -  |iL  .  p,  vhere  p  is  a  constant  of  integxation. 
r 

If  k  x  *  i,  then 


C  x  v  s  C  1  v  /, 


v  =  A -  ju  £  +  p|  ( Hodograph) 


Let  us  show  that 


p  «  u  e  i  . 


-  P  =  C  xv+|i£s(rxv)xv+nl, 
r  r 


p  =  (v  .  v)  r  •  (r  •  v)  V  -  |1  £ 

r 


(n) 
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If  p  is  constant,  it  may  be  evaluated  at  any  point  in  the  orbit, 
and  we  will  choose  the  pericenter.  Then 


P  =  rp  i  -  n  £  =  (v^  rp  -  u)  3*  =  U  e  £ 


as  we  had  to  prove, 
Using  (ll) 


ff  s  (V  .  $)  ?  - 


(?  .  ?)  ?  -  u  I 
r 


and  applying  Poisson's  formula,  it  is  easily  found  that 
p  =  2(v  •  If)  r  -  (r  •  R)  v  -  (r  •  v)  R 
Now,  in  an  equatorial  system 

/  cos  w  cos  ft  +  sin  u  sin  ft  cos  I 
i  =  I  cos  w  sin  ft  +  sin  w  cos  ft  cos  I 
\  sin  u  sin  I 

Therefore 

p  =  -£(uer)=u|eT4*ei^ifi+e£!iw+e2iil 
dt  *“  dft  dw  d* 

From  p  .  i  =  G  u  e  ,we  obtain  (10)  in  a  new  form 


e 


(12) 


(13) 


and  from  (12) 

»  2(7  .  if)  <?  •  n  -  (r  •  (f)  (?  •  t)  -  (?  •  7)  <R  •  T)  (14) 

From 

p  •  ^  *  ue  (cos  I  ft  -f  w) 
dw 


we  deduce 


w  =  -  cos  I  ft  + 


-  hi 

llE. 

ue 


(15) 
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Therefore 


a  =  -  -.?£S.,I  h  X  &  1 ,»  j.  [ 2g  .fyh  •£)+ 
2C  sin2  I  ^  on^f16^  >>  du  / 


In  resume,  equations  (6),  (8),  (9),  (10)  or  (14),  (16)  give  the 
variational  equations  for  a,  n,  I,  e,  u. 

4.  Time  of  Pericenter 


-1/2  3/2  ,  v 

T  =  t-  |i  a  (E-e  sin  E) 

-Applying  Poisson's  formula 

t  =  -  n1/s  grad-  £a3/2  (E-e  sin  E)  !•  R  (17) 

Now 

gradr*  a  =  — ^  v 

grad-  E  = - i - \lP™?  ("(hr2  +  C2)  v  -  h(r  •  v)  rl  -  SL  v  - 

e  sin  E*»(i2eL  JpJ 


grad^  e  =  «JL.  (hr2  +  C2)  v  -  h(r  •  ?)? 

H2e  L  J 

and  these  relations  enable  us  to  compute  t,  from  (17).  Alternatively 
we  may  compute  a  where  a  =  l  -  nt. 

SUMMARY 

The  equations  for  the  variation  of  constants  have  been  deduced  in 
a  more  direct  way,  without  the  necessity  of  computing  Lagrange 
brackets.  On  the  other  hand  it  seems  that  the  present  method 
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is  suitable  for  any  integral  and  gives  the  direct  variation  of  a 
desired  element  by  the  use  of  a  single  equation.  Furtheremore,  it 
gives  the  variations  of  the  integrals  of  motion  in  an  elegant 
way.  It  is  hoped  that  the  preeent  method  simplifies  the  deduction 
tc  a  great  extent. 

The  method  is  essentially  equivalent  to  a  method  suggested  to  the 
author  by  Dr.  C.  J.  Cohen,  where  Poisson's  variational  equation  is 
applied  directly  to  the  elliptic  elements,  considered  as  functions 
of  position  and  velocity. 
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